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Fermi has discovered two giant gamma-ray-emitting bubbles that extend nearly 10 kpc in diam- 
eter. We propose that periodic star capture processes by the galactic supermassive black hole, Sgr 
A*, with a capture rate < 10 -5 yr _1 and energy release ~ 10 52 erg per one capture can produce 
shocks in the halo, which accelerate electrons to the energy 1 TeV. These electrons generate radio 
emission via synchrotron radiation, and gamma-rays via inverse Compton scattering with the relic 
and the galactic soft photons. Estimates of the diffusion coefficient from the observed gamma-ray 
flux explains consistently the necessary maximum energy of electrons and sharp edges of the bubble. 



I. INTRODUCTION 

The recently discovered Fermi bubbles are symmet- 
ric gamma-ray structures derived from the Fermi LAT 
data in the energy range 1-100 GeV. The bubbles elon- 
gate above and below the Galactic plane for about 8 
kpc and their radius is about 3 kpc . Observations 
show a very sharp outer boundary of the bubble. The 
gamma-ray intensity sharply drops outward the bub- 
bles Q. 

The origin of the bubble is still enigmatic and up 
to now a few models were presented in the literature. 
Our group assumed that the Fermi bubbles originated 
from star capture events which occurred in the GC ev- 
ery 10 4 — 10 5 years @. These events form giant shocks 
propagating through the central part of the Galactic 
halo and thus produce accelerated electrons with en- 
ergies < 10 TeV whose scattering on background pho- 
tons is responsible for the bubble gamma-ray emission. 

Processes of particle acceleration by the bubble 
shocks in terms of sizes of the envelope, maximum 
energy of accelerated particles, etc. may differ sig- 
nificantly from those obtained for SNs that may lead 
to the maximum energy of accelerated protons much 
larger than can be reached in SNRs. In this respect, 
we assume that acceleration of protons in Fermi bub- 
bles may contribute to the total flux of the Galactic 
cosmic rays (CR) above the 'knee' break (> 10 15 eV). 



II. BUBBLE HYDRODYNAMICS 

We assume in Q that the central black hole cap- 
tures a star every t ~ 10 4 — 10 5 years. As a result 
the total energy £ ~ 10 52 erg releases in the Galac- 
tic center in the form of 100 MeV proton which heat 
the central 20 pc up to the temperature ~ 10 keV. 
This heating produces a shock propagating into the 
surrounding medium. In the simplest case this situa- 



tion can be described by a solution obtained by [3] for 
the adiabatic explosion in the exponential atmosphere 
with the density profile p(z), 



(1) 



p(z) = p exp | 



where z is the coordinate perpendicular to the Galac- 
tic plane. For parameters of the Galactic halo po = 
0.25 cur 3 and zq = 1 kpc. The shock propagating 
in to the halo forms in the exponential atmosphere a 
double bubble structure elongated in z-direction. The 
radius of the bubble at the height z and at the time t 
is 



r = 2zo arccos 
where 
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(3) 



V is a current volume bounded by the shock 



a(t) 



V(t) = 2tt J r 2 (z,t)dz, 
o 

a is the position of the shock top 



a(t) = -2z In ( 1 - ^ ) . 



(4) 



(5) 



7 is the polytropic coefficient, and a and A are num- 
bers. 

For the finite time t\ determined from the condition 
y{t\) = 2zo the shock breaks through the exponential 
atmosphere and the bubble top a(ti) tends to infinity 
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while the bubble radius in the Galactic plane (z = 0) 
tends asymptotically to the the value 



2zn arccos (1/2) ~ 2 kpc , 



(6) 



that is comparable with the radius of Fermi bubbles. 
For So ~ 10 52 erg and po ~ 0.25 cm -3 the value of t\ is 
about 3 x 10 8 yr. We want to remark that t\ sensitively 
depends on the injected energy and the density profile 
of the halo. For example if the injected energy is £q ~ 
3 x 10 52 erg and po ~ 0.1 cm -3 , t\ can be reduced by 
nearly an order of magnitude. 

Then for a periodic star capture the bubble interior 
is filled with shocks propagating in series one after 
another through the bubble interior stopping at the 
radius ~ 2z$ That gives formally a stationary sideway 
boundary. 

The realistic situation has to be described by a 
set of dissipative hydrodynamic equations, which take 
into account the shocks propagation in non-uniform 
medium and various dissipation processes including 
shock heating, energy transfer into cosmic rays, slow- 
ing down due to accumulating material etc. These 
processes are ignored in the the Kompaneetz solution. 
Shocks should disappear when their speed is lower the 
local sound speed, shocks should disappear when their 
speed is lower the local sound speed. Then the sideway 
boundary of the Bubble is simply given by T\, ~ v s tdi S 
where v s is the sound speed and tdis is the charac- 
teristic time of the shock dissipation because of e.g. 
particle acceleration at the shock front. 



III. ELECTRON ACCELERATION AND 
GAMMA-RAYS FROM THE BUBBLE 

We assume that the bubble gamma-rays are pro- 
duced by IC scattering of electrons on the relic pho- 
tons. For the rate of synchrotron and inverse Comp- 
ton energy losses dE/dt — f3E 2 the maximum energy 
of electrons accelerated by shocks estimated in 

the Bohm diffusion limit is 



El 



leHu 2 
3cf3 



(7) 



that gives e.g. for the shock velocity u 



10 8 



s , the magnetic field strength H = 10 G and the 
energy density of relic photons w p h = 0.25 eV cm -3 
the maximum energy of accelerated electrons about 
E e max ~ 5 x 10 13 eV. 

In Fig. [T] we show the expected spatial distribu- 
tions of gamma-ray emission from the bubble for sin- 
gle shock and multiple shocks cases. From this figure 
one can see that the single shock model is unable to 
reproduce the data. However, for the parameters of 
star capture model these data are nicely described. 

The expected spectrum of gamma-ray emission 
from the Bubble due to IC scattering of the electrons 
is shown in Fig. [2] 
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FIG. 1: Spatial distribution of the gamma-ray emission. 
Data are from Q. Top: in case of single shock. Dotted 
line correspond to D=10 29 cm 2 /s, solid to D=10 30 cm 2 /s, 
dashed to D=10 31 cm 2 /s and dash-dotted D=10 32 cm 2 /s, 
Bottom: in case of several shocks distributed in accordance 
with ©. 
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FIG. 2: The spectrum of gamma-ray emission from Fermi 
bubble in case of multi-shock acceleration. Data points 
are taken from []]|. 



IV. PROTON ACCELERATION IN THE 
BUBBLE AND THE ORIGIN OF THE "KNEE" 
COSMIC RAYS 

CRs within with energies below E ~ 10 15 eV are 
generally attributed to SNRs in our Galaxy. We as- 
sume that some of the CRs produced by SNRs in the 
Galactic disk are re-accelerated by shocks in the Bub- 
ble to energy above 10 15 eV that explains the origin 
of CRs beyond the knee. 

For the multi-shock structure in the bubble an av- 
erage distance L between separate shocks given by 



L = tqu = 30 



tq 



s 3 x 10 4 yr / \ 10 8 cm/s 
where u is the shock front speed. 



pc. 



(8) 
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FIG. 3: CR spectrum derived in the model of particle re- 
acceleration in the bubble. The data are summarized in 



5]. 



If the value of L exceeds the scale of particle accel- 
eration by a single shock which is Id ~ D/u where 
D is the spatial diffusion coefficient near a shock and 
u is the shock velocity, then particle acceleration by 
shocks is pure stochastic, which describes by a mo- 
mentum diffusion coefficient (see Q) 



and particle momentum, and Q(p,z,p) describes CR 
injection by supernova remnants in the disk with en- 
ergies E < 10 15 eV with the spectrum Q oc p~ 4 . The 
spectrum of CRs injected by SNRs and re-accelerated 
in the bubble is shown in Fig. [3] 



V. CONCLUSION 

We have shown that series of shocks produced by 
a sequential stellar captures by the central black hole 
can further re-accelerate the protons emitted by SNRs 
up to energies above 10 15 eU. The predicted CR spec- 
trum contributed by the Bubble may be E~ u where 
v ~ 3 for I0 15 eV < E < 10 19 eV that explains the 
knee CR spectrum. 

The regime of electron acceleration in the bubble 
is quite different from that of protons. It is a com- 
bination of single and multishock accelerations. In 
this case we have a cut-off of the electron spectrum 
at E > 3 10 13 eV and flattening of the spectrum at 
E < 100 GeV that explains nicely the bubble gamma- 
ray spectrum and the sharp edge spatial distribution 
observed by [l[ if this emission is due to IC on the 
relic photons. 
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Then the equation describing particle production by 
SNRs in the disk, their re-acceleration in the bubble 
and propagation in the Galaxy can be presented in 
the form 



0_ 

dz 
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'p2Q^ \*(J>,P)P Yp ' =-Q(P' Z 'P)> 



(10) 



where p and z are the cylindrical spatial coordinates, 
p is the particle momentum. D(p,p) is the spatial dif- 
fusion coefficient, which is a function of coordinates 
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